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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE j 

Group Art Unit 21 54 j 
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| 



DECLARATION UNDER 37 C.F.R. § 1.131 

Commissioner for Patents 

United States Patent and Trademark Office 

P.O. Box 1450 

Alexandria, VA 22313-1450 

Dear Sir: 

1. This declaration is to establish the status of the invention in the above- 
captioned U.S. Patent Application prior to March 26, 1993, the filing date of U.S. Patent 
No. 5,368,026, issued on November 29, 1994, to Swedlow et al. ("Swedlow"). 

2. We understand that Claims 1-10, 12-16, 18-25, 27-36, and 38-46 of the 
above-captioned U.S. Patent Application Ser. No. 10/677,050 (the '050 Patent 
Application), which was filed September 30, 2003, titled "Signal Processing Apparatus," 
have been rejected in the February 27, 2008 Office Action as being unpatentable over 
Swedlow. 

3. We are inventors of the subject matter of the above-referenced claims in 
the '050 Patent Application. 

4. At least prior to March 26, 1993, we conceived of and reduced to practice 
a signal processing apparatus such as the device and method described and claimed in 
the above-referenced claims of the '050 patent application. 

5. As evidence thereof, we provide herewith Exhibit A which is a design 
note ("the Design Note") describing the subject matter of the above-referenced claims. 
The Design Note was created and is dated prior to March 26, 1993. We conceived the 
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ideas referenced therein prior to the creation of the Design Note. Any references to 
dates on the Design Note have been redacted. 

6. At Masimo Corporation ("Masimo"), the term "Stork" refers generally to 
Masimo's pulse oximetry technologies. 

7. At least prior to March 26, 1993, we compiled, loaded onto hardware, and 
used a software application called "STORK" for its intended purpose, that is, to measure 
blood oxygen saturation. 

8. As evidence thereof, we provide herewith Exhibits B and C which are 
screen captures of a directory of files, including ".c" source code files, ".h" header code 
files, and ".gdd" debugging code files compiled into the application file STORK. As 
shown on the screen captures, each of these files includes a last edited date prior to 
March 26, 1993, which has been redacted herein. 

9. Exhibits D and E are the header and source code files, titled "vector.h," 
and "vector.c," respectively, written in the C programming language, showing the 
mathematical functions for calculating a vector cross correlation. 

10. Exhibits F and G are the header and debugging files, titled "stork.h," and 
"stork.gdd," respectively, written in the C programming language. The header file 
defines the variable for cross correlation and the debugging file displays the value of 
cross correlation for instructive purposes. 

11. Exhibit H is the source code file, titled "iawa.c," written in the C 
programming language, which calculated the cross correlations. 

12. Exhibit I is the source code file, titled "satfil.c," written in the C 
programming language, which qualifies bins of signal data responsive to intensity 
signals from a detector monitoring a patient, the qualification including determining 
whether the calculated cross correlation is within expected limits. 

13. Exhibit J is the source code file, titled "ui.c," written in the C programming 
language, which is the user interface for the STORK application and which indicates an 
expected limit of cross correlation. 

14. The foregoing files and other exhibits demonstrate that we constructed 
and tested an embodiment and performed processes meeting every element of the 
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above-referenced claims prior to March 26, 1993 sufficient to demonstrate that the j 

constructed embodiment operated for its intended purpose. J 

15. We hereby declare that acts leading to the reduction to practice of the j 
subject matter claimed in the '050 Patent Application were performed in the United 

States. j 

I declare that all statements made herein of my own knowledge are true and that j 
all statements made on information and belief are believed to be true, and further that 

these statements were made with the knowledge that willful, false statements and the j 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 

Title 1 8 of the United States Code and that such willful, false statements may jeopardize j 

the validity of the application or any patent resulting therefrom. ] 



Dated: By: 



I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



Mohamed K. Diab 



Dated: 



By: 



Walter M. Weber 
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I declare that all statements made herein of my own knowledge are true and that 
all statements made on information and belief are believed to be true, and further that 
these statements were made with the knowledge that willful, false statements and the 
like so made are punishable by fine or imprisonment, or both, under Section 1001 of 
Title 18 of the United States Code and that such willful, false statements may jeopardize 
the validity of the application or any patent resulting therefrom. 



Dated: 



By: _Al 




Massi E. Kiani 
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The goal of this note is to propose an artifact detector for the Stork algorithm. Specifically, the artifact 
detector is designed such that it assigns a real number to each snapshot of data that it processes, and then 
compares that real number to a predetermined threshold value. Snapshots of data which yield artifact 
detector values that lie below the threshold value are said to contain significant amounts of artifact. 
Snapshots of data (hat yield artifact detector values which lie above the threshold value are said to contain 
insignificant amounts of artifact. 

In principle, the Stork algorithm accepts as input two signals which are of the form (after log and 
bandpass) 

S a (t) = s a (t) + n a (t) 

(1) 

S b (0 = s b (t) + n b (t) 

where 

s b (t) - m x 83(1) 

(2) 

i«b(0-«>2 «•»«>. 

Any input signals which are not expressable in the above form are said to consist of artifacts. 
Consequently to determine whether or not a signal is artifact free, it is sufficient to determine whether or 
not equations (1) and (2) are satisfied. 

^----~T5neway to check whether or not equations (1) and (2) have been met is to 
C__ between S a (t) and S b (t). Specifically one could compute the cross correlation 



f S a (t) S b (t) dt = f < s a (t) s b (t) + n a (t) n b (t) ) dt + j ( s a (t) n b (t) + s^t) n a (l) ) dt 



which assuming that the signals, s a (t) and s b (t), are uncorrected with the noises, n a (t) and n b (t), may be 
expressed as 



J S a (t) S b (t) dt = J ( s a (t) s b (t) + n a (t) n b (t) ) dt. 



A normalized cross correlation may be computed as follows 

J S a (t) S b (t) dt I ( s a (t) s b (t) + n a (t) n b (t) ) dt 

[/ S2 a (t) dt J S2 b(t ) dt]0-5 (J ( s 2 a(t) + „2 a{t) } dt f ( ^(t) + „2 b(t) } dtJ 0.5' 



The solutions which give the extrema of this cross correlation functional are given by equation (2). Any 
input signals which are not expressible in the form given in equation (2) yield reduced normalized 
correlation values. In Figure I, these extrema of the normalized correlation functional are plotted vs. the 
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"Red" Signal - to - Noise ratio and co,/co 9 . One should note that the normalized correlation is not 
necessarily unity and precaution must be taken, particularly when the difference between arterial 
saturation and the background saturation is large. 

One may rigorously prove that the magnitude of the normalized cross correlation is always less than or 
equal to unity by use of the inner product 

x-y = fx(t)y(t)dt. 
Using the above notation, the normalized cross correlation may be written as 

V% + n a*"b 
[(s a .s a + n a .n a )(s b .s b + n b »ntflp-S 

or equivalently 

__ i a l s i 

tS a -S a S b .S b j°-5 ' 

But, according to the Schwarz Inequality 

(S a .S b )2 Z (S a .S a )(S b .S b ). 
Hence the normalized correlation is always between 1 and -1. 
ALGORITHM IMPLEMENTATION 

1. Compute the quantity 

2S a S b 

"Es2 a z S 2 b] oT" 

2. Compare quantity to a threshold, say 0.75 

THINGS THAT STILL NEED TQ BE INVESTIGATED 

Investigate performance vs. snapshot size 
Investigate performance vs signal/noise correlation 
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4KB csourcefte 
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4KB CSourcefte 

4KB csourcefte 

2KB CHeaderfie 

13KB csourcefte 

SKB CSourcefte 

1KB CHeoderrte 

4KB CHeada-fle 

2KB CHeaderfie 

4KB csourcefte 

4KB CSourcefte 

7KB Fife 

3KB csourcefte 



3KB CHeaderfie 
9 KB GOD Hie 
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14KB csourcefte 
10 KB CSourcefte 
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File: vector. h 

Description: 

Collection of routines that work on vectors. 
Masimo Standard Library. 

Public Functions: 
min_index 
max_index 
innerjoroduct 
remove_mean 



Returns index to minimum element in vector 
Returns index to maximum element in vector 
Returns point by point product of vectors 
Removes the mean value from each point in a vector 

Computes Dot product of two vectors 
Computes Correlation Coefficient 



Private Functions: 



01/15/93 
01/25/93 
02/25/93 
02/25/93 



RJM 
RJM 
RJM 
RJM 



Initial Version 
add remove_mean 
add VECTOR_DotProduct 
Add VECTOR_CrossCore0 



#include <masimo.h> 



extern INT32 
min_index ( 

FLOAT32 

INT32 



extern INT32 
max_index { 

FLOAT32 

INT32 

extern VOID 
inner_product ( 
FLOAT 3 2 
FLOAT 3 2 
INT32 
FLOAT32 

extern VOID 
remove_mean ( 
FLOAT32 
FLOAT32 
INT32 



*data, 
size) ; 



*data, 
size) ; 



*vl, 
*v2, 
size, 
*v3) ; 



*in, 
*out, 
size) ; 



extern FLOAT32 
VECTOR_DotProduct ( 
FLOAT32 *u, 
FLOAT32 *v, 
INT32 size) ; 

extern FLOAT32 
VECTOR_CrossCorr0 ( 
FLOAT32 *u, 
FLOAT 3 2 *v, 
INT32 size) ; 



/* input, pointer to data 

/* input, number of data points 



/* input, pointer to data 

/* input, number of data points 



input , 
input , 
input , 



first vector 
second vector 
size of vectors 



input , data 

output, data with mean removed 
input, number of points 



input, first vector 
input, second vector 
input, size of vectors 



input, first vector 
input, second vector 
input, size of vectors u < 
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File: vector. c 
Description: 

Collection of routines that work on vectors. Designed for 
Masimo Standard Library. 

Public Functions : 

max_index 
inner_jjroduct 
remove_mean 
VECTOR_Dot Produc t 
VECTOR_CrossCorrO 



Returns index to smallest element in vector 
Returns index to largest element in vector 
Computes the inner product of two vectors 
Removes the mean value from each point in a vector 
Compute Dot Product of two vectors 
computes cross correlation time=0 



Private Functions : 



01/15/93 RJM Initial Version 
02/25/93 RJM Add VECTOR_DotProduct 
02/25/93 RJM Add VECTOR_CrossCorrO 



#define MODULE_ID 

#include <masimo.h> 
# include <math.h> 
#include "vector. h" 



/* platform definitions 
/* sqrt function 
/* self 



INT32 

min_index ( 
FLOAT32 
INT32 

INT32 
INT32 
FLOAT32 



*data, 
size) { 



j; 

tmp; 



/* input, pointer to data 

/* input, number of data points 



BUG1 (data) ,- BUG1 (size) ; 



j = 0; 

tmp = data [j] ; 

for (i = l; i < size; i++) { 
if (data[i] < tmp) { 
tmp = data [i] ; 



} 



j - i; 



} 

INT32 

max_index ( 
FLOAT32 
INT32 

INT32 
INT32 
FL0AT32 



*data, 
size) { 



j; 

tmp; 



/* input, pointer to data 

/* input, number of data points 



BUG1 (data) ; BUG1 (size) ; 
j = 0; 

tmp = data[j] ; 

for (i = l; i < size; i++) { 
if (data[i] > tmp) { 
tmp = data [i] ; 

return ( j ) ; 



VOID 

innerjproduct ( 
FLOAT32 
FLOAT32 
INT32 
FLOAT32 

INT32 



*vl, 
*v2, 
size, 
*v3) { 



/* input, first vector 

/* input, second vector 

/* input, size of vectors 

/* output, v3 = vl * v2 



for (i = 0; i < size; i++) { 
v3 [i] = vl[i] * v2 [i] ; 

} 



VOID 

remove_mean ( 
FLOAT 3 2 
FLOAT 3 2 
INT32 



*out, 
size) { 



input , data 

output, data with mean removed 
input, number of points 



BUG1 (in) ; BUG1 (out) ; BUG1 (size) ; 
mean = 0.0; 

for (i = 0; i < size; i++) mean += in[i]; 
mean /= (FLOAT32) size; 

for (i = 0; i < size; i++) out [i] = in[i] 



/* sum all points */ 



FLOAT 3 2 

VECTOR_Dot Product ( 
FLOAT32 *u, 
FLOAT32 *v, 
INT32 size) { 

FLOAT32 sum; 



/* input, first vector 
/* input, second vector 
/* input, size of vectors 



BUG1 (u) ; BUG1 (v) ; BUG1 (size) ; 
sum = 0.0; 

for (i = 0; i < size; i++) { 
sum += u [i] * v [i] ; 

} 

return (sum) ; 



135 } 



FL0AT32 

VECTOR_CrossCorrO ( 
FLOAT32 *u, 
FLOAT 3 2 *v, 
INT32 siz 

FLO AT 3 2 den; 



0 { 



/* input, first vector 
/* input, second vector 
/* input, size of vectors u and v 



BUG1 (u) ; BUG1 (v) ; BUG1 (size) ; 

den = SQRTF (VFCTORJDotProduct (u, u, size) * 
VECTORJDotProduct (v, v, size)); 

if (den ==0.0) { 

C = 1.0; 
} else { 

c = VECTOR_DotProduct (u, v, size) / den; 
return (c) ; 
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File: stork. h 

Description: Stork header file included by ALL stork modules 
Public Functions: 
Private Functions: 



History: 




RJM 


11/06/92 


RJM 


12/08/92 


RJM 


12/16/92 


RJM 


01/04/93 


RJM 


01/05/93 


RJM 


01/14/93 


RJM 


01/15/93 


RJM 


01/16/93 


RJM 


01/18/93 


RJM 


01/18/93 


RJM 


01/19/93 


RJM 


01/25/93 


RJM 


01/25/93 



Initial version (vl.29) 
Prevent multiple inclusion of stork. h 
Change vsinc params to all floats 
Add MAX_BINS 

Add seed to bin structure 

Put conditional code around "C" only stuff 
Make memory pool sizes platform dependant 
Change sat result structure for filtering 
Make HR a vector 
Add more sat params 

Add hr_max_delta and sat_max_delta params 
Make SAT equation adjustable 
Make HR algorithms switchable 



#define VERSION 1.541 
#define VSINC_LEAN 0 
#define ADC_SNAPSHOT_SIZE 



/* floating point software version number 
/* Generate code for Lean Architecture 



Memory Allocation 



#if (PLATFORM . 



1) 



#define MEM_STACK_SIZE 
#define MEM_POOL_SIZE 
#define DEBUG_MEM_S I ZE 
#define GDD_BUF_SIZE 
#define DBG_FIFO_SIZE 

#else 

#define MEM_STACK_SIZE 
ttdefine MEM_POOL_SIZE 
#define DEBUG_MEM_SIZE 
ftdefine GDD_BUF_SIZE 
#define DBG_FIFO SIZE 



32768 
512 

20000 
100 



26000 
512 
6000 
100 



/* data memory 

/* data memory 

/* data memory 

/* pmda 

/* pmda 



[was 6000) 



#endif 

/* The following contant controlls ALL debug code for the STORK 
/* beta software. No required code should be conditioned on this 
/* constant. 



ttdefine DEBUG 1 



/* enables debug code 



/* 



The following constants control what data is captured in the dbg_fifo. 
The dbg_fifo is controlled by the ADC interrupt routine. 



#define DEBUG_SIM 
#define DEBUG_ADC 
#define DEBUG_LED 
#define DEBUG_SUMS 
ttdefine DEBUG_RAW 
#define DEBUG HSL 



/* simulate ADC data */ 

/* stores raw adc conversions in fifo */ 

/* stores tick mark for LED change */ 

/* stores the reading for each state */ 

/* stores raw ir and raw red */ 

/* stores High Speed Link data */ 



/* 

Low Speed Link Control 
*/ 

#define LSL_CHECK_OVERRUN 1 /* check HC11 FIFO */ 

#define LSL_GDD_GROUP_SIZE 6 /* number of gdd packets per lsl_tick */ 

/* __. 

High Speed Link Control 
*/ 

#define HSL_RECORD 0 /* records ADC data after ambient subtract */ 

#define HSL_PLAYBACK 0 /* uses data from HSL instead of ADC */ 

#define HSL_TEST_PATTERN 0 /* generate rotating l's pattern */ 

#define IMPULSE 0 /* generate impulses for usrl */ 

/* -- -- 

The following constants control the ANC. 

*/ 

#define ANC_FIXED_CELL 0 /* fixed cell normalized anc */ 

/* 

System Wide Type definitions 

*/ 

#define FIXED_CELL 1.0 /* anc_types */ 

#define ACCORD IAN 2.0 /* anc_types */ 

/* The following constants describe the size of the following */ 

/* structures. */ 

#define MAX_BINS 6 
#define NUM_HR_ALG 7 

#define BIN_RESULTS_SIZE 8 

#define SAT_PARAMS_SIZE 22 

#define SAT_RESULTS_SIZE 14 

#define HR_PARAMS_SIZE 7 

#define HR_RESULTS_SIZE 8 

ttifndef ASSEMBLY 

typedef struct { 

INT32 num_points; /* number of input points */ 

FL0AT32 anc_type; /* type of anc */ 

FL0AT32 lambda; /* nANC parameter */ 

FLOAT32 min_error; /* nANC parameter */ 

INT32 anc_settle_time; /* nANC parameter */ 



1 


FLOAT32 


num_cells ; 






FLOAT 3 2 


peak limitl 






FLOAT 3 2 


peak_epsilon; 


140 


FLOAT32 


*bp_coef ; 






INT32 


bp size; 






FLOAT32 


search_f irst ; 




FLOAT32 


search step; 




INT32 


search_size 




145 


INT32 


bin_count ; 






FLOAT32 


cov limit; 






FLOAT3 2 


DVHS init; 






FLOAT 3 2 


DVHS0 init; 






FLOAT 3 2 


DVAST timer 




150 


FLOAT 3 2 


DVAST delta 






FLOAT32 


DWST timer 






FLOAT32 


DWST delta 






FLOAT3 2 


max sat delta; 




INT32 


sat eqn; 




155 


FLOAT32 


dual_lambda ; 




FLOAT32 


correlation; 




} SAT_PARAMS ; 








typedef struct 


{ 




160 


FLOAT32 


ir_rms ; 






FLOAT32 


rd_rms ; 






FLOAT32 


cov; 






FLOAT32 


arterial sat; 




FLOAT32 


venous_sat ; 




165 


FLOAT32 


seed_sat ; 






FLOAT 3 2 


nuinjpeaks ; 






FLOAT32 


correlation 






} BINJRESULTS; 






170 


typedef struct 


{ 






FLOAT32 


Spa02_f inal 






FLOAT32 


Spv02_f inal 






FLOAT32 


Spa02 unfil 






FLOAT32 


Spa02_fill; 




175 


FLO AT 3 2 


Spa02 dual; 






FLO AT 3 2 


DVHS; 






FLOAT 3 2 


DVHS0 ; 






FLOAT3 2 


DVAST; 






FLOAT32 


DVAST time; 




180 


FLOAT 3 2 


DWST; 






FLOAT 3 2 


DWST time; 






FLOAT 3 2 


art_omega_pr ; 




FLOAT32 


ven_omega_pr ; 




FLOAT32 


time; 




185 


BIN RESULTS 


bins [MAX BINS] ; 




} SAT RESULTS; 







typedef struct { 

190 FLOAT 3 2 resolution; 

FLOAT 3 2 raw_sampl ing_rate ; 

FLOAT 3 2 min_hr; 

FLOAT32 max_hr; 

FLOAT 3 2 nst; 

195 FLOAT 3 2 max_hr_delta; 

INT32 hr_enable; 

FLOAT 3 2 hr_lambda; 
} HRJPARAMS; 

200 typedef struct { 

FLOAT 3 2 hr [NUM_HR_ALG] ; 

FLOAT 3 2 hr_final; 

FLOAT 3 2 BPM_per_bin; 

INT3 2 num_bins ; 





nANC parameter 


*/ 


/* 


peak picker parameter 


*/ 


/* 


peak picker parameters 


*/ 


/* 


pointer to filter coef . 


*/ 


/* 


number of filter coef . 


*/ 


/* 


first saturation value . 


*/ 


/* 


saturation step size. 


*/ 


/* 


number of steps 


*/ 




number of bins 


*/ 


/* 


coyariance limit 


*/ 


/* 


initial value for DVHS 


*/ 


'. 


initial value for first DVHS 


*/ 


1 * 


seconds 


*/ 


/* 


change for DVAST 


*/ 


/* 


seconds 


*/ 


/* 


change in DWST 


*/ 


/* 


max change in sat w/motion 


*/ 


/* 


saturation equation number 


*/ 


/* 


dual ANCs as VSINC1 inputs 


*/ 


/* 


bins OK only if >= corr. 


*/ 



/* 


signal rms 


*/ 


/* 


signal rms 


*/ 


/* 


covariance 


*/ 


/* 


resulting arterial sat 


*/ 


/* 


resulting venous sat 


*/ 


/* 


seed saturation 


*/ 


/* 


peaks in power curve 


*/ 


/* 


raw dataset correlation 


*/ 



/* resulting arterial sat 

/* resulting venous sat 

/* unfiltered arterial sat 

/* first stage of filtering 

/* dual ANCs as input to VSINC1 

/* default value high sat 

/* default value high sat first 

/* default value art sat (PR) 

/* time DVAST was last updated 

/* default value ven sat (PR) 

/* time DWST was last updated 

/* arterial omega for pulse rate 

/* venous omega for pulse rate 

/* time for these results 

/* bin results 



/* HR resolution in BPM */ 
/* raw data sampling rate in Hz */ 
/* minimum hr allowed (BPM) */ 
/* maximum hr allowed (BPM) */ 
/* noise_spec thresshold */ 
/* maximum change in HR w/motion*/ 
/* bits control with alg */ 
/* ANC algo. uses this lambda */ 

/* (BPM) from each alg */ 

/* (BPM) after all filtering */ 

/* BPM per bin of spectrum */ 

/* number of bins per spectrum */ 



205 ^} 



Exhibit G 



; Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 
; 

5 ; File: stork. gdd 







Description: Debug Display Control File 


10 




Public Functions: 








Private Functions: 








Notes : 








Channels 0 to 39 are defined in SPT.H 


15 




Channels 97 to 88 are for status line 






History: 








RJM 12/10/92 


Initial version 


20 




RJM 12/21/92 Add algorithm parameter 






RJM 12/22/92 Add bin results table and plots 






RJM 12/28/92 Add new HR algorithm results j 






RJM 12/30/92 


Change to one algorithm 




; 


RJM 12/3 0/92 Add peak_epsilon parmaters 


25 




RJM 01/05/93 Add Seed sat to bin data, remove others 






RJM 01/16/93 New SAT results structure 






RJM 01/18/93 Add dvhs init and dvhsO init parameters 






RJM 01/18/93 Add hr5 result 






RJM 01/18/93 Add more sat params 


30 




RJM 02/22/93 Real-Time data Plotting 






RJM 02/23/93 Clean up Channel Assignments 






Debugl screen 


i 




p 


3,2 


display key = 1 




DEBUG SCREEN 1 






30 


050, 025, 0375, 125, 7 


raw data 




34 


050,175,0375,275,5 


analysis data 


40 


33 


050,325,0375,425,5 


plethy 




32 


050,475, 0375,575, 5 


freq spect j 




20 


775,025,1000,100,5 


power curve 




21 


775,150,1000,225, 0 


bin 1-5 ; 




22 


775,230,1000,305, 0 




45 


23 


775,310,1000,385,0 






24 


775,390,1000,465,0 






25 


775,470,1000,545,0 






31 


775,600,1000,725,5 


ADC Snapshot ! 




00 


400,005,0000,000,0 


system params j 


50 


01 


575, 005, 0000, 000, 0 


Sat Results 




03 


575,225, 0000, 000, 0 


HR Results 




04 


575, 370, 0000, 000, 0 


Calib Results j 




02 


010,610,0000,000,0 


bin results j 




17 


400,490,0014,000,12 


ADC OOR count j 


55 


98 


000,741,0014,000,12 


'Working' Text I 




99, 000,741, 0015, 000, 12 


Version* | 




97, 000,755, 0015, 000, 12 
END 


file name. . 


60 




Debug2 screen Raw, Spectrum, Plethy, tables 




P 1,1 


display key = 2 




DEBUG SCREEN 2 






30 


050,025,425, 150,5 


raw 


65 


32 


050,225,425,330,7 


freq 




33 


050,400,425,525,5 


plethy 




00 


450,005,000,000,0 


System Parameters 




01 


625, 005, 000, 000, 0 


Sat Results 



03,625,225,0000,000,0 
04,625,400,0000,000,0 
98,050,570,014,000,12 
99,050,570,015,000,12 
97,050,585,015,000,12 



HR Results 
Calib Results 
'Working' Text 
Version# 
file name. . 



Debug3 screen power curve, bin graphs 



P 2,1 

DEBUG SCREE 
20, 050,020, 
21, 050, 180, 
22, 050,330, 
23,450,020, 
24,450,180, 
25,450,330, 
02, 060,460, 
98,050,590, 
99,050,590, 
97,350,590, 
END 



380,130,7 
380,280,7 
380,430,7 
780, 130, 7 
780,280,7 
780,430,7 
000, 000, 0 
014,000,12 
015,000,12 
015, 000, 12 



bin results 
' Working ' Text 
Versions 
file name. . 



Debug4 screen Raw, Spectrum, Plethy, tables 



P 3,0 
DEBUG i 
30,050, 
32,050, 
33,050, 
00,300, 
01,470, 
03,470, 
04,470, 
98, 000, 
99, 000, 
97, 000, 



270,125,5 
270,280,5 
270,430,4 
000,000,0 
000,000,0 
0000, 000, 0 
0000, 000, 0 
014,000,12 
015,000,12 
015,000,12 



,- display key = 4 

raw 
freq 
plethy 

system params 
sat results 
HR Results 
Calib Results 
' Working ' Text 
Version* 
file name . . 



Debug5 screen power curve, bin graphs 



P 4, 0 
DEBUG 
20, 050 
21, 050 
22,050 
23,375 
24,375 
25,375 
02,020 
98, 020 
99, 020 
97,250 
END 



SCREEN 5 
020,295, 
135,295, 
250,295, 
020,620, 
135,620, 
250,620, 
345,000, 
471, 014, 
471, 015, 
471,015, 



100,7 
215, 7 
320,7 
100,7 
215,7 
320,7 
000, 0 
000, 12 
000, 12 
000,12 



display key = 5 



bin results 
' Working ' Text 
Version* 
file name. . 



Customer Screen #1 

P 5,0 

Customer Screen 

52,000,000,015,000,28 

77,050,050,450,150,21 

78, 050,200,450,300,21 

53, 500, 080, 014, 000, 16 

54,500,230,014,000,16 

15,500,110,014,000,16 

16,500,260,014,000,16 

31,500,310,630,400,00 

17, 500,415, 011, 000, 08 



display key = 6 

stork legend 

Real time IR Data (23) 

Real time RD Data 

Spa02 : 

Pulse : 

Spa02 value 

Pulse value 

ADC Snapshot 

ADC OOR count HGK 



18,500,430,011,000,08 
14,500,445,011,000,08 
40,500,460,012,000,08 



Correlation 
Tick 

calibrate: 1/0 message 



Customer Screen #2 

P 6,0 

Customer Screen 

52,000,000,015,000,28 

35,050,050,450,150,21 

33,050,200,450,300,05 

53,500,080,014,000,16 

54,500,230,014,000,16 

15,500,110,014,000,16 

16,500,260,014,000,16 

31,500,310,630,400,00 

17,500,415,011,000,08 

18,500,430,011,000,08 

14,500,445,011,000,08 

40,500,460,012,000,08 



stork legend 

raw analysis portion 

plethy 

Spa02 : 

Pulse : 

Spa02 value 

Pulse value 

ADC Snapshot 

ADC OOR count HGK 

Correlation 

Tick 

calibrate: 1/0 message 



Front End Screen 

P 7,1 

Front End Screen 

04, 500, 010, 0000, 000, 0 

31,500,200,775,575,7 

17,050,050,014,000,12 

77,050,200,450,375,23 

78,050,400,450,575,23 

40,050,020,011,000,28 



Calib Results 
ADC Snapshot 
ADC OOR count 
Real time IR Data 
Real time RD Data 
calibrate: 1/0 message 



Channel Definitions 



Usrl snapshot 



Raw Data 
Point number 
Amplitude 



G 31 

ADC SNAPSHOT 
Point Number 
Amplitude 



Frequency Spectrum 



G 32 

Spectrum 
HR (BPM) 
Amplitude 



Master Power Curve 



G 20 

Power Curve 
Sat 



G 21 
Bin #1 
Sat 
Amp 



G 22 
Bin #2 
Sat 
Amp 



Bin #3 

Sat 

Amp 



G 24 
Bin #4 
Sat 



G 25 
Bin #5 
Sat 

240 Amp 

; Plethysmograph 

G 33 

Clean Plethysmograph 
245 Time 
Amp 

; analysis data 

G 34 



; raw analysis data 

255 G 35 

Raw Data (before Masimo) 

Index 

Amp 



260 



Info Strings 



F 98,0 

Working 

F 99,0 

265 Stork Version %6.3f 



Algorithm Parameters 



T 0,1,32 
270 ALG 1 

cov_limit,%9.1f 
num_points , %9 . Of 



anc_type, %9 . Of 
lambda, %9.4f 
min_error, %9.4f 
anc_settle, %9 . Of 
num_cells,%9.0f 
peak_liraitl, %9.2f 
peak_epsilon, %9 . 2f 
search_f irst , %9 . 2f 
search_step, %9 . If 
search_size , %9 . Of 
bin_count , %9 . Of 
sampling rate, % 9. Of 
dvhs_init, %9. If 
dvhs0_init,%9.1f 
DVAST_timer,%9.0f 
DVAST_delta,%9.1f 
DWST_timer,%9.0f 
DWST_delta, %9 . If 
max_sat_delta, %9 . If 
max_hr_delta, %9 . If 
sat egn,%9.0f 
hr enable, %9. Of 
Cal Target, % 9. 4 f 
Cal Err Band,%9.4f 
Cal Set. Time, %9. Of 
Cal IR PGA,%9.2f 
Cal RD PGA,%9.2f 
hr lambda, %9.4f 
dual lambda, % 9. 4 f 
Correlation, %9.4f 

Sat Results 

T 1,1,14 
ALG 1 

>Spa02_f inal , %9 . 4f 
>Spv02_f inal , %9 . 4f 
>Spa02_unf il , %9 . 4f 
>Spa02_fill,%9.4f 
>Spa02_dual,%9.4f 
DVHS, %9.4f 
DVHSO, %9.4f 
DVAST,%9.4f 
DVAST time,%9.4f 
DWST, %9.4f 
DWST time,%9.4f 
>art_omega_pr , %9 . 4f 
>ven_omega_pr , %9 . 4 f 
>time, %9.4f 

HR Results 

T 3,1,8 
Alg 1 

>hrO (ir),%9.4f 
>hrl (sig_ref ) , %9.4f 
>hr2 (deriv) ,%9.4f 
>hr3 (ir/n),%9.4f 
>hr4 (sig/n) , %9 .4f 
>hr5 (spec5) ,%9.4f 
>hr6 (ANC),%9.4f 
>hr_final,%9.4f 

Bin Results 

T 2,6,8 
Master Bin 
Bin 1 



Bin 2 
Bin 3 
Bin 4 
Bin 5 

ir_rms,%9.4f 
rd_ntis, %9 .4f 
>cov, %9.4f 
>arterial sat,%9.2f 
>venous sat,%9.2f 
>seed sat,%9.2f 
>num_jpeaks , %9 . Of 
>corr,%9.6f 

Calib Results 

T 4,1,10 
C Results 
ir_trans, %9 .4f 
rd_trans,%9.4f 
desired_ir,%9.4f mA 
desired_rd,%9.4f mA 
actual_ir,%9.4f mA 
actual_rd,%9.4f mA 
ir__gain, %9.2f 
rd_gain, %9.2f 
ir_delta,%9.4f 
rd_delta,%9.4f 

; display strings 
; 

M 40,0 
Calibrating 
S 52,1 

Masimo Stork 

S 53,0 

Spa02 : 

S 54,0 

Pulse : 

F 14,0 

Tick %12.0f 

F 15,0 

%3.0f %% 

F 16,0 

%3.0f BPM 

F 17,0 

ADC OOR Count %8.0f 
F 18,0 

Sensor Pos. %8.4f 

; Real-Time Plotting 

R 77,1,400 

Plethysmographies Data (before Masimo) 
Time 

Amplitude 
R 78,1,400 

Plethysmographies Data (after Masimo) 
Time 

Amplitude 
X 12, 77, 0, 9 
IR 

X 13,78,0,12 



end of items 



Exhibit H 



Copyright (c) Masimo Corporation (tm) 1992, 1993. 

File: iawa.c 

Description: 

Industry Algorithm Without Averaging 

Public Functions: 
init_iawa ( ) ; 
iawa_sat ( ) ; 

Private Functions: 

Notes : 



All Rights Reserved. 



RJM 11/05/92 Based on s.c {VI. 29) 

RJM 11/20/92 use param and result structures 

RJM 12/08/92 Make compatible with VSINC 

RJM 01/05/93 Compute separate seed for each bin 



#def ine MODULE_ID 


1015 




# include 


<masimo.h> 


/* 


platform specification 


*/ 


# include 


<memory.h> 


/* 


special memory allocation 


*/ 


#include 


< debug. h> 


/* 


assert , grab 


*/ 


# include 


<math.h> 






# include 


" stork. h" 


/* 


system parameters 


*/ 


# include 


"fir.h" 


/* 


block fir module 


*/ 


# include 


"sateqn.h" 


/* 


saturation equations 


*/ 


# include 


"iawa.h" 


/* 


self 


*/ 


# include 


"vector .h" 


/* 


cross corr 


*/ 


# include 


"gddcli.h" 


/* 


General Debug Display Client 




# include 


"spt .h" 







#if (PLATFORM == 1) 
extern VOID 
presentation ( 

FL0AT32 *data, 
INT32 first, 
INT32 num, 
INT32 code) ; 

#endif 
VOID 

init_iawa (void) { 

^ /* nothing yet */ 

VOID 

compute_log ( 

FLOAT 3 2 *in, 
FLOAT32 *out, 
INT32 n) { 



/* this BS is to allow UI.C code to store. 



/* input data 
/* output data 
/* number of points 



/* This routine computes the natural log of the input signal */ 

/* Because of the filtering in the demodulator, there are times */ 

/* when the resulting signal can go below zero, this log routines */ 

/* limits small values so that the log is gracefully limited. */ 

#define LOG00001 -11.51293 



INT32 



BUG1 (in) ; BUG1 (out) ; BUG1 (n) ; 



for (i = 0; i < n; i++) { 

if (in[i] > 0.00001) out [i] = LOGF (in [i] ) ; 
else out[i] = LOG00001; 

} 

} 

VOID 

iawa_sat ( 

FLOAT 3 2 *log_ir, /* input, normalized ir data */ 

FLOAT32 *log_rd, /* input, normalized rd data */ 

SAT_PARAMS *sp, /* input, parameters */ 

SAT_RESULTS *sr) { /* output, results */ 



INT32 


b; 


/* 


bin counter 


*/ 


INT32 


bin_size ; 


/* 


point in a bin 


*/ 


INT32 


f_size; 


/* 


size of filtered data 


*/ 


BIN RESULTS 


*r; 


/* 


points to results for each bin 


*/ 


FLOAT32 


*f ir; 


/* 


filtered ir 


*/ 


FLOAT32 


*f_rd; 


/* 


filtered rd 


*/ 


FLOAT32 


* xaxis; 


/* 


for green marker lines 


*/ 



BUGl(log_ir) ; BUGl(log_rd) ; BUGl(sp); BUGl(sr); 

assert (log_ir != NULL) ; 
assert (log_rd != NULL) ; 
assert(sp != NULL) ; 
assert(sr != NULL) ; 

f_size = sp->numjpoints - sp->bp_size; 

bin_size = (f_size - sp->anc_settle_time) / sp->bin_count; 
assert (f_size > 0) ; 

mark_mem ( ) ; 

f_ir = get_mem(f_size) ; 
f_rd = get_mem(f_size) ; 
xaxis = get_mem(l000) ; 

/* band-pass filter IR and RD signals */ 

blkfir (log_ir, sp->num_points, sp->bp_coef, sp->bp_size, f_ir) ; 

blkfir (log_rd, sp- >num_points , sp->bp_coef, sp->bp_size, f_rd) ; 

/* output analysis data */ 

b=0; /* build a dataset that demarks the bins divisions */ 

assert (f_size-sp->anc_settle_time < 1000); 
for(i=0; i < f_size - sp->anc_settle_time; i++) 

xaxis [i] = (FLOAT32) (b-5000) ; 
if ( ! (i%bin_size) ) b=10000-b; 



j = sp->anc_settle_time- l; 
b = f_size - j ; 

GDDCLI_Chart(SPT_AnalysisData,3,b, 0.0, 1.0); 
GDDCLI_Data (SPT_AnalysisData, b, &f_rd [ j ] ) ; 
GDDCLI_Data (SPT_AnalysisData,b, &f_ir [j] ) ; 
GDDCLI_Data (SPT_AnalysisData,b, xaxis) ; 

#if (PLATFORM == 1) 

presentation ( f _rd, j ,b, 0) ; 
#endif 



calculate correlations */ 

sr->bins [0] .correlation = VECTOR_CrossCorrO ( &f _ir [ j ] , &f_rd[jj , f_size-j) ; 
for (b = 1; b <= sp->bin_count; b++) 
r = &(sr->bins [b] ) ; 

r->correlation = VECTOR_CrossCorrO (&f_ir [j] , &f_rd [j ] ,bin_size) 
j+= bin_size; 

Square both signals */ 

for (i = 0; i < f_size; i++) 

f_ir[i] *= f_ir[i]' ; 
f_rdti] *= f_rd[i] ; 

Integrate using binning */ 

sr->bins [0] . ir_rms = 0.0; 
sr->bins [0] . rd_rms = 0.0; 

sp->anc_settle_time - 1; 

for (b = 1; b <= sp->bin_count; b++) 

r = &(sr->bins [b] ) ; 

r->ir_rms = 0.0; 
r->rd_rms = 0.0; 

for (i = 0; i < bin size; i++) 
{ 

r->ir_rms += f_ir[j]; 

r->rd_rms += f_rd[j] ; 

) d "" 

/* update master bin */ 

sr->bins [0] .ir_rms += r->ir_rms; 
sr->bins [0] .rd_rms += r->rd_rms; 

r->ir_rms = SQRTF (r- >ir_rms / bin_size) ; 
r->rd_rms = SQRTF (r- >rd_rms / bin_size) ; 

r->seed_sat = calc_sat (r->rd_rras / r->ir_rms, sp->sat_eqn) ; 



/* compute master bin results */ 
r = &(sr->bins [0] ) ; 

r->ir_rms = SQRTF (r- >ir_rms / f_size) ; 
r->rd_rms = SQRTF (r- >rd_rms / f_size) ; 
r->seed_sat = calc_sat (r->rd_rms / r->ir_rms, sp->sat_eqn) ; 
free_mem() ; 



/* was bin_size HGK */ 



Exhibit I 



Copyright (c) Masimo Corporation (tm) 1992. All Rights Reserved. 
File: satfil.c 

Description: Spa02 output filtering algorithms 



Public Functions: 
reset_satf il 
sat_filter 

Private Functions: 
max_dif f_alg 



Resets all filters 
runs all sat filters 



ROM 10/09/92 Initial version 

RJM 10/12/92 Allow multiple filtering streams 

ROM 11/06/92 Bring up to coding standards 

RJM 01/16/93 Merge binning algorithms, DVHS filters, etc 

RJM 01/18/93 If COV is low, and 1 peak don't report Sv02 

RJM 01/26/93 Limit final Spa02 by highest seed of all bins 

RJM 01/26/93 Compare max-diff values to Spa02_final 



ftdefine MODULE_ID 1004 



ttinclude <masimo.h> 
#include " stork. h" 
#include "satfil.h" 
#include "sateqn.h" 
# include <memory.h> 
#include <debug.h> 
#include <vector.h> 
#include <math.h> 

VOID 

reset_satf il ( 

SAT_PARAMS *sp, 
SAT_RESULTS *sr) { 

BUGl(sp) ; BUGl(sr) ; 

assert (sp != NULL); 
assert (sr != NULL) ; 



sr->DVHS = sp->DVHS_init; 
sr->DVHS0 = sp->DVHS0_init; 
sr- >Spa02^f inal = 0.0; 
sr- >Spa02_unf il = 0.0; 
sr->Spa02_f ill = 0.0; 
sr->Spv02_f inal = 0.0; 
sr- >Spa02_dual = 0.0; 
sr->DVAST = 0.0; 
sr->DWST = 0.0; 
sr->DVAST time : 



/* platform description file */ 

/* system parameters */ 

/* self * */ 

/* calc_omega */ 

/* Special memory routines */ 

/* assert, grab */ 

/* vector functions */ 

/* */ 



/* input, saturation parameters 
/* inout, saturation results 



sr->DWST tit 



■sp->DVAST_timer • 
= -sp->DWST_timer • 



65 /* 



Name: 



max_dif f_alg 



This routine finds the two bins with the largest difference 
between arterial and venous sat and of these two returns 
the saturation which is closest to the previous sat. For 
more detail see Stork design note 27. 



1 . Select the bins that have a good arterial and venous 
saturation and order these bins based on the difference between 
the venous and the arterial saturation. 

2. Select the two bins with the greatest difference. If there 
are no bins resulting from step 1. 

3 . Of these two bins select the one closest to the previous 
Spa02_final value. 

4. If there is no previous value, then return the largest value. 



static FLOAT 3 2 
max_diff_alg( 



"sat PARAMS 


*sp, 


/* 


input, saturation parameters */ 


SATJRESULTS 


*sr) { 


/* 


inout, saturation results */ 


INT32 


i; 






BIN RESULTS 


*b; 


/* 


general bin ptr */ 


INT32 


*good bins; 


/* 


array of good bin numbers */ 


FLOAT32 


*diffs; 


/* 


array of differences */ 


INT32 


num good bins; 


/* 


number of good bins */ 


FLOAT32 


sat; 


/* 


return value */ 


FLOAT32 


satl; 


/* 


arterial sat of bin with largest diff*/ 


FLOAT32 


sat2 ; 


/* 


sat of bin with 2nd largest diff. */ 



BUG1 (sp) ; BUGl(sr) ; 



mark_mem ( ) ; 

good_bins = get_mem (sp- >bin_count) ; 
diffs = get_mem(sp->bin_count) ; 

/* select bins that have good arterial and venous sat */ 
/* R154 : bin must also have good correlation HGK */ 

num_good_bins = 0; 

for (i = 1; i <= sp->bin_count; i++) { 
b = &(sr->bins [i] ) ; 

if (b->arterial_sat > 0 && b- >venous_sat > 0 && 
b->correlation >= sp- >correlation) 



{ 



good_bins [num_good_bins] = i; 

dif f s [num_good_bins] = b->arterial_sat 

num_good_bins++ ; 



b- >venous_sat ; 



/* Find Two bins with largest difference */ 

switch (num_good__bins) { 

case 0: 

sat = sr->Spa02_unf il; 
break; 



case 1: 



i = max_index (dif f s, num_good_bins) ; 

sat = sr->bins [good_bins [i] ] . arterial_sat ; 

break; 

default : 

/* find bin with largest diff */ 

i = max_index (dif f s, num_good_bins) ; 

satl = sr->bins [good_bins [i] ] .arterial_sat; 

/* find bin with second largest diff */ 

diffs[i] = 0.0; 

i = max_index(diffs ( num_good_bins) ; 

sat2 = sr->bins [good_bins [i] ] .arterial_sat; 

/* If there is no previous Spa02 choose highest sat */ 

if (sr->Spa02_f inal ==0.0) { 
sat = satl; 

if (sat2 > satl) sat = sat2; 
} else { 

/* Otherwise, choose one closest to Spa02_final */ 
sat = satl; 

if (FABSF (satl - sr->Spa02_final) > 

FABSF ( sat2 - sr->Spa02_final) ) sat = sat2; 



f ree_mem ( ) ; 
return ( sat ) ; 



sat_f ilter 

This routine takes the master power curve results and the 
five bins and computes an overall arterial and venous 
saturation based on the following algorithm. For more 
detail see Stork design note 27. 



/* max sat for final filter output 
/* min sat for final filter output 



r) { 



_f ilter ( 
SAT_PARAMS 
S AT_RE SULTS • 



BIN_RESULTS *mb; 
FLOAT32 delta; 
INT32 i; 
FLOAT 3 2 max seed; 



master bin ptr 

maximum seed from all bins 



BUGl(sp) ; BUGl(sr) ; 



asserttsp != NULL); 



assert (sr != NULL); 

/* If master bin cov is less than limit then use its sat value */ 
/* otherwise, use the sat of the max_diff_alg of the bins. */ 

mb = &(sr->bins [0] ) ; 

if (mb->cov < sp->cov_limit) { 

sr->Spa02_unf il = mb->arterial_sat; 

/* Update DVHS */ 

if (mb->arterial_sat > sr->DVHS0) { 
sr->DVHS = mb->arterial_sat; 
sr->DVHS0 = sr->DVHS; 

} 

/* DVAST algorithm */ 

if (mb->arterial_sat > 0.0 && mb->arterial_sat < sr->DYHS) { 
sr->DVAST = mb->arterial_sat; 
sr->DVAST time = sr->time; 

} 

} else { 

sr->Spa02_unf il = max_dif f_alg (sp, sr) ; 
/* DWST algorithm */ 

if (mb->venous_sat > 0.0 && mb->num_jpeaks > 1) { 
sr->DWST = mb->venous_sat; 
sr->DWST time = sr->time; 

, 1 

/* first stage of filtering (DVHS) */ 
sr->Spa02_fill = sr->Spa02_unf il; 

if (sr->DVHS < sr->Spa02_fill) sr- >Spa02__f ill = sr->DVHS; 

/* second stage of filtering */ 
/* If there is motion and a previous value, limit the change */ 
/* to max_sat_delta . */ 

if (mb->cov > sp->cov_limit && sr- >Spa02_f inal > 0.0) { 

delta = sr->Spa02_f ill - sr- >Spa02_f inal ; 

if (delta > sp- >max_sat_delta) delta = sp- >max_sat_delta; 

if (delta < - sp- >max_sat_delta) delta = - sp- >max_sat_delta; 

sr->Spa02_f inal += delta; 
} else { 

sr->Spa02_f inal = sr->Spa02_f ill; 

/* Find largest seed from all bins */ 

max_seed = mb->seed_sat; 

for (i = 1; i <= sp->bin_count; i++) { 

^ if (sr->bins [i] .seed_sat > max_seed) max_seed = sr->bins [i] . seed_sat; 

if (sr->Spa02_f inal < max_seed) sr->Spa02_f inal = max_seed; 
if (sr->Spa02_final > MAX_SPA02) sr->Spa02_f inal = MAX_SPA02; 
if (sr->Spa02_final < MIN_SPA02) sr->Spa02_f inal = MIN_SPA02; 

/* Venous logic */ 

if (mb->cov < sp->cov_limit && mb- >num_peaks < 2) { 
sr->Spv02_f inal =0.0; 



} else { 

sy->Spv02_f inal = mb->venous_sat; 

/* Compute arterial and venous sat for use by Pulse Rate algorithms */ 
sr- >art_omega_pr = 0.0; 

if (sr->time - sr->DVAST_time < sp->DVAST_timer) { 

if ( FABSF ( sr- >DVAST - sr- >Spa02_f ill) < sp->DVAST_delta) { 
sr->DVAST time = sr->time; 

} 

sr- >art_omega_pr = calc_omega (sr->DVAST, sp->sat_eqn) ; 
} else { 

if (sr->Spa02_fill > 0.0) { 

^ sr->art_omega__pr = calc_omega (sr->Spa02_f ill, sp->sat_eqn) ; 
sr- >ven_omega_pr = 0.0; 

if (sr->time - sr- >DWST_time < sp->DWST_timer) { 

if (FABSF (sr->DWST - mb- >venous_sat) < sp- >DWST_delta) { 
sr->DWST time = sr->time; 

} 

sr->ven_omega_pr = calc_omega (sr->DWST, sp->sat_eqn) ; 
} else { 

if (mb->venous_sat > 0.0 && 

mb->cov > sp->cov_limit && mb- >num_peaks > 1) { 
sr->ven_omega_pr = calc_omega (mb->venous_sat, sp->sat_eqn) ; 

} 

} /* satfil */ 



310 #ifdef STANDALONE 

ftinclude <stdlib.h> 
#include <stdio.h> 

315 VOID main (void) { 



} 

320 #endif 



Exhibit K 
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I 



/* 

Copyright (c) Masimo Corporation (tm) 1992, 1993. All Rights Reserved. 



File: ui.c 
Description: 

This is the STORK software user interface module. This platform 
independant user interface provides the ability to start and stop 
the oximeter, change system parameters, and control the input data. 
There are minor differences in these features listed below: 



Feature 



PC Platform 



DSP Platform 



Start/Stop 
Change system param 
LED calibration 
PAN through data 
Rerun same data 
Standard Mode 

Public Functions: 
init_ui 

user_interface 

Private Functions: 
store_raw 
store_param 
store_results 
single_step 
test ^pattern 

Notes : 





History: 




35 


RJM 


12/09/92 




RJM 


12/28/92 




RJM 


12/30/92 




RJM 


12/30/92 




RJM 


01/18/93 


40 


RJM 


01/18/93 




RJM 


01/18/93 




RJM 


01/25/93 




RJM 


01/25/93 




RJM 


01/25/93 


45 


RJM 


02/02/93 




RJM 


02/11/93 




RJM 


02/12/93 




RJM 


02/19/93 




RJM 


02/23/93 


50 


RJM 


02/24/93 




HGK 


03/08/93 



Yes 
Yes 
Yes 



Yes 
Yes 
Yes 



Standard initialization 
Start the user interface 



send raw data via LSL 

send parameters via LSL 

send results via LSL 

Run a single snapshot of data 

send a test pattern via LSL 



Initial Version based on main.c 

New HR module vsi_hr 

New default values for parameters 

Add peak_epsilon parameter 

Add dvhs_init and dvhs0_init parameters 

Add hrfil module 

Keep track of time 

Set default for DVHS0 to 95 

Make SAT egn parameter 

hr_enable param 

Add command for DSP Platform version 

Validate data after get_snapshot 

New default values for coefs 

Integrate file list feature from HK. 

Add Calibration parameters 

Improve ADC Snapshot 

move data_j?tr code to usrl_c 



ttdefine I 


!ODULE_ID 


1024 




#include 


<masimo .h> 


/* 


Platform description 


*/ 


#include 


<pckeybrd . h> 


/* 


key definitions 


*/ 


#include 


<memory . h> 


/* 


memory manager 


*/ 


# include 


<debug.h> 


/* 


virtual platform debugger 


*/ 


# include 


<math . h> 








# include 


"sateqn.h" 


/* 


saturation equations 


*/ 


# include 


"stork.h" 


/* 


System parameters 


*/ 


# include 


"adc.h" 


/* 


Demodulation module 


*/ 


# include 


"usrl .h" 


/* 


Filtering and Decimation 


*/ 


#include 


"calib.h" 


/* 


LED calibration 


*/ 


#include 


"vsinc .h" 


/* 


VSI noise cancellation module 


*/ 


#include 


"satfil.h" 


/* 


Saturation Filters 


*/ 


#include 


"hrfil.h" 


/* 


Pulse Rate Filters 


*/ 





„ . , , 


"lsl h" 


/* 


low s eed link module 


* / 


70 


#incl de 


" ' awa h" 


/* 


I d st 1 ' th ' th t 






# include 


"post h" 




post band^as^f ilter ° U avera 9 in 9 






# include 


"vsihr .h" 


/* 


VSI Heart Rate Module 


*/ 




# include 


"vsipl .h" 


/* 


VSINC Clean photoplethysmograph 


*/ 




# include 


" gddcli .h" 


/* 


General Debug Display Client 


*/ 


75 


#include 


11 spt .h" 


/* 


Stork Packet Types 


*/ 




# include 


"vsidual.h" 


/* 


dual ANC sat version 


*/ 




# include 


"ui.h" 


/* 


Self 


*/ 




#if (PLATFORM == 1) 









#include <stdio.h> 
static BOOL presentation_f lag; 
ttendif 



typedef struct { 

FLOAT3 2 min_val; 

FLOAT32 max_val; 

FLOAT3 2 inc_val ; 

FLOAT32 def_val; 
} PARAM_DESC; 



/* minimum value 

/* amount to change 
/* default value 



typedef struct { 





SAT_PARAMS 


*vp; 


/* 


pointer to vsinc params 


*/ 


I 




SAT RESULTS 


*vr; 


/* 


pointer to vsinc results 


*/ 


\ 


95 


FLOAT32 


*power ; 


/* 


pointer to power curve 


*/ 






HR PARAMS 


*vhp; 


/* 


pointer to Heart Rate params 


*/ 






HR RESULTS 


*vhr ; 


/* 


pointer to Heart Rate results 


*/ 






FLOAT32 


*f spec; 


/* 


frequency spectrum 


*/ 






FLOAT32 


*pleth; 


/* 


clean plethysmograph 






100 


CALIB PARAMS 


*cp; 


/* 


pointer to calibration params 


*/ 






CALIB RESULTS 


*cr; 


/* 


pointer to calibration results 


*/ 






} ALGORITHM SPEC; 












105 


/* The following constants 


control the user adjustable */ 




\ 




/* parameters. 






*/ 




if 




#define MAX ALGORITHMS 




1 








110 


#define NUM_PARAMS 




32 






| 



#define 
#define 
#define 
#def ine 
#define 
#def ine 
#define 
#def ine 
#define 
#define 
#define 
#define 
#define 
#def ine 
#def ine 
#define 
#define 
#define 
#define 
#define 
#def ine 
#define 
ftdefine 
#define 
#define 



P_cov_limit 

P_num_points 

P_anc_type 

P_lambda 

P_min_error 

P_anc_settle_timi 

P_num_cells 

P_peak_limitl 

P_peak_epsilon 

P_search_f irst 

P_search_step 

P_search_size 

P_bin_count 

P_sampling_rate 

P_dvhs_init 

P_dvhs0_init 

P_DVAST_timer 

P_DVAST_delta 

P_DWST_timer 

P_DWST_delta 

P_max_sat_delta 

P_max_hr_delta 

P_sat_egn 

P_hr_enable 

P_cp_target 



ftdefine P_cp_error_band 
ttdefine P_cp_settle_time 
#define P_cp_ir_gain 
#define P_cp_rd_gain 
#define P_hr_lambda 
#define P_dual_lambda 
ttdefine P_correlation 

static PARAMJDESC params [NUM_PARAMS] = { 



/* 








min _ 




max 




anc 


def */ 


/* 


cov_limit 


0 


* / 


o 




200 


' 




0, 


40 . 0 


/* 


num_points 


1 


* / 


100 




1000 




10 


0, 


650.0 


/* 


anc_type 


2 














0, 


1.0 


/* 


lambda 


3 


*/ 








n 




0001 


0.75 


/* 


min_error 


4 


* / 












0001 


0.0005 


/* 


anc settle time 5 


* / 


Q 




1000 






o, 


149.0 


/* 


num cells 


6 


* / 






50 






o , 


12. 0 


/* 


peak_limitl 


7 


* / 


o 










01, 


0. 02 


/* 


peak_epsilon 


8 


* / 


o 










01, 


0. 03 


/* 


search_f irst 


9 






0 ' 


120 






25, 


104 . 0 


/* 


search step 




* / 












10, 


0.5 


/* 


search_size 


11 


* / 






300 






o. 


195.0 


/* 


bin_count 


12 


* / 












o , 


5 . 0 


/* 


sampling rate 


13 


* / 


10 


0 ' 


200 


' 




o, 


62.5 


/* 


dvhs init 


14 


* / 


40 


0 ' 


105 




0 


1, 


99.0 


/* 


dvhsO_init 


15 


* / 


40 




105 




0 


1, 


95.0 


/* 


DVAST timer 


16 




Q 


0 ' 


100 


Q ' 


1 


o, 




/* 


DVAST delta 


17 


*/ 


0 


o', 


100 


0, 


0 


1, 


2.0 


/* 


DWST timer 


18 


*/ 


0 


o, 


100 


o, 


1 


0, 


10.0 


/* 


DWST_delta 


19 


*/ 


0 


o, 


100 


o, 


0 


1, 


4.0 


/* 


max_sat_delta 


20 


*/ 


0 


o, 


20 


0, 


0 


1, 


2.0 


/* 


max_hr_delta 


21 


*/ 


0 


o, 


40 


0, 


0 


1, 


4.0 


/* 


sat_eqn 


22 


*/ 


0 


0, 


1 


o, 


1 


o, 


1.0 


/* 


hr_enable 


23 


*/ 


0 


o, 


255 


o, 


1 


0, 


251.0 


/* 


cp_target 


24 


*/ 


0 


o, 


1 


o, 


0 


05, 


0.25 


/* 


cp_error_band 


25 


*/ 


0 


o, 


0 


5, 


0 


01, 


0.10 


/* 


cp_settle_time 


26 


*/ 


0 


o, 


300 


o, 


2 


00, 


75.0 


/* 


cp_ir_gain 


27 


*/ 


0 


0, 


4096 


o, 


0 


1, 


2.0 


/* 


cp_rd_gain 


28 


*/ 


0 


o, 


4096 


o, 


0 


1, 


2.0 


/* 


hr_lambda 


29 


*/ 


0 


1, 


1 


o. 


0.0001, 


0.9999 


/* 


dual_lambda 


30 


*/ 


0 


1, 


1 


o. 


0.0001, 


0.9999 


/* 


correlation 


31 


*/ 


0 


1, 


1 


0, 


0.0001, 


0.9999 



param_val [MAX_ALGORITHMS] [NUM_P ARAMS ] ; 



static FLOAT32 
VOID 

init_ui (void) { 

/* nothing yet */ 

store raw,plethy data to ASCII files to prepare for presentations 
#if (PLATFORM == 1) 
VOID 

presentation ( 
FLOAT32 
INT 3 2 
INT32 
INT32 



{ 



*data, 
first, 



code) 



static CHAR 



filename [65] ; 
str [65] ; 
*fp; 
i; 



if ( !presentation_f lag) return; 
switch (code) 



. 1 

#endif 



printf ("Enter B/A filename (.B,.A will be added) 
scanf ("%s", filename) ; 
sprintf (str, "%s.B", filename) ; 
goto part2; 

sprintf (str, "%s . A" , filename) ; 

if ( (fp=fopen(str, n w+") )==NOLL) break; 
for (i=0;i<num;i++) { 

f printf (fp, "%f\n» ,data [f irst+i] ) ; 

fclose(fp) ; 
break; 



static VOID 

store_raw ( 

SAT_PARAMS 
FLOAT32 
FLOAT 3 2 
INT32 
INT32 



size, 

time_stamp) { 



/* input, parameters */ 

/* input, ir datav */ 

/* input, red data */ 

/* input, number of points */ 

/* input, starting point number */ 



BUG1 (sp) ; BUG1 (ir) ; BUG1 (rd) ; BUG1 (size) ; BUG1 (time_stamp) ; 
#if ( PLATFORM == 1) 



markjnem ( ) ; 

t = get_mem(size) ; 

assert(t != NULL); 

j = sp->bp_size/2 + sp- >anc_settle_time; 

k = size - sp->bp_size - sp->anc_settle_time; 



/* 1st pt */ 
/* #usable pts */ 



for(i=0;i<j ;i++) t[i] = ■ 
for( ,-i<j+k ;i++) t[i] = 
for( ;i<size;i++) t [i] = 
GDDCLI_Chart (SPT_UsrlSnapShot, 3, size, ( FL0AT3 2 ) t ime_s tamp , 1.0); 
GDDCLI_Data(SPT_UsrlSnapShot, size, ir) ,- 
GDDCLI_Data(SPT_UsrlSnapShot, size, rd) ; 
GDDCLI_Data(SPT_UsrlSnapShot, size, t ); 



GDDCLI_Chart ( SPT_RawAnalysisData , 
GDDCLI_Data (SPT_RawAnalysisData, 
GDDCLI_Data (SPT_RawAnalysisData, 



2,k, (FLOAT32)time_stamp, 1.0); 
k, &ir[j]) ; 
k, &rd[j]) ; 



static VOID 
store ^parameters ( 

INT32 num_algorithms) { 

INT32 i; 

BUG1 (num_algorithms) ; 

for (i = 0; i < num_algorithms; i++) { 



GDDCLIJData (SPT_SystemParam, NUM_PARAMS, param_val [i] ) ; 



static VOID 
store_results ( 

ALGORITHM_SPEC 

INT32 

SAT_PARAMS 
SAT_RESULTS 
HR_PARAMS 
HR_RESULTS 
FLOAT32 
FLOAT32 
FLOAT32 
INT32 
INT32 
FLOAT32 



*alg, 

num_algorithms) { 

*vp; 
*vr; 
*vhp ; 
*vhr ; 



*power ; 
*f spec; 
*pleth; 
i; 

j ; 



Frequency SPECtrum 
clean pleth. 



BUGX (alg) ; BUG1 (num_algorithms) ; 
mentry("store_results") ; 
/* frequency spectrum */ 
i = 0; 

vp = alg[i] .vp; 
vr = alg [i] ,vr; 
power = alg [i] .power; 
vhp = alg[i] .vhp; 
vhr = alg[i] .vhr; 
fspec = alg[i] .fspec; 
pleth = alg[i] .pleth; 

/* frequency spectrum */ 

GDDCLI_Chart(SPT_FreqSpec, 1, vhr->num_bins, vhp->min_hr, 

vhr->BPM_j>er_bin) ; 
GDDCLI_Data(SPT_FreqSpec, vhr->num_bins, alg [i] . fspec) ; 

/* power curve */ 

GDDCLI_Chart (SPT_MasterPowerCurve, num_algorithms, vp->search_size, 

vp->search_first, - (vp->search_step) ) ; 

for (i = 0; i < num_algorithms; i++) { 

GDD CL I_Da t a ( S PT_Ma s t e r Powe r Curve , vp->search_size, alg [i] .power) ; 

for (j =0; j < vp->bin_count; j++) { 

GDDCLI_Chart{SPT_PowerCurvel + j, num_algorithms, vp->search_size, 

vp->search_first, - (vp->search_step) ) ; 

for (i = 0; i < num_algorithms ; i++) { 

GDDCLI_Data (SPT_PowerCurvel + j, vp->search_size, 

&(alg[i] .power [(j+1) *vp->search_size] ) ) ; 



/* plethysmographic signal */ 

GDDCLI_Chart (SPT_CleanPleth, num_algorithms , pleth_size (vp) , 0.0, 1.0); 
for (i = 0; i < num_algorithms; i++) { 

GDDCLI_Data(SPT_CleanPleth, pleth_size (vp) , alg [i] .pleth) ; 



#if (PLATFORM == 1) 

presentation (alg[0] .pleth, 0 , pleth_size (vp) ,1) ; 
#endif 



/* store tabular results */ 

for (i = 0; i < num_algorithms; i++) { 

GDDCLI_Data(SPT_SatResults, SAT_RESULTS_SIZE, (FLOAT32 * ) alg [i] . vr) ; 
^ GDDCLI_Data(SPT_HRResults, HR_RESULTS_SIZE, (FLOAT32 *> alg [i] .vhr) ; 

for (i = 0; i <= alg[0] .vp->bin_count; i++) { 

GDDCLIJData (SPTJBinResults, BIN_RESULTS_SIZE, 
^ (FLOAT32 *)&(alg[0] .vr->bins[i] )) ; 

/* master bins' correlation sent separately (times 100) */ 

temp = alg [0] .vr->bins [0] .correlation * 100.0; 
GDDCLI_Data (SPT_correlation, 1, &temp) ; 

mexit("store results"); 

} 

static INT32 
compute_snapshot_size ( 

ALGORITHM_SPEC *alg, 

INT32 num_algorithms) { 

INT32 i; 
INT32 size; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

/* find largest snapsnot size */ 
size = 0; 

for (i = 0; i < num_algorithms; i++) { 

if (alg [i] . vp->num_j>oints > size) size = alg [i] . vp- >num_j>oints ; 

return (size) ,- 

} 

static VOID 
new_jpatient ( 

ALGORITHM_SPEC *alg, 

INT32 num_algorithms) { 

INT32 i; 

BUG1 (alg) ; BUG1 (num_algorithms) ; 

for (i = 0; i < num_algorithms ; i++) { 
reset_satf il (alg [i] .vp, alg [i] .vr) ; 
reset_hrfil(alg[i] .vhp, alg[i] .vhr) ; 



static VOID 
single_step ( 

ALGORITHM_SPEC *alg, /* input, algorithm structure */ 

INT32 num_algorithms, 

INT32 raw_snapshot_size) { 

INT32 i ; /* algorithm */ 

FLOAT32 *ir; /* raw ir data */ 





FLOAT 3 2 


*rd; 


/* 


raw rd data. 


*/ i 


410 


FLOAT32 


*log ir; 


/* 




*/ 




FLOAT32 




/* 


log of ird data 


*/ 1 




FLOAT32 


*adc _ snap- 
_ 


/* 


snapshot of adc data 


*/ 




SAT P ARAMS 


*vp; 






1 


415 


SAT RESULTS 


*vr; 






j 




HR PARAMS 


*vhp; 






\ 




HR RESULTS 


*vhr; 






I 




FLOAT32 


♦power; 










FLOAT32 


*fspec; 


/* 


Frequency SPECtrum 


*/ 


420 


FLOAT32 


*pleth; 


/* 


clean pleth. 


*/ 




INT32 


time; 









BUGl (alg) ; BUG1 (num_algorithms) ; 



mentry ( " single_step" ) ; 
markjnem ( ) ; 



BUGl (raw_snapshot_size) , 



/* Allocate data required for current algorithms 

ir = get_mem (raw_snapshot_size) 

rd = get_mem (raw_snapshot_size) 

log_ir = get_mem(raw_snapshot_size) 

log_rd = get_mem(raw_snapshot_size) 

adc_snap = get_mem (ADC_SNAPSHOT_SIZE) ; 

ADC_GetSnapshot (adc_snap) ; 
adc_snap [0] = -1.0; 
adc_snap[l] = 1.0; 

GDDCLI_Chart(SPT_ADCSnapShot, 1, ADC_SNAPSH0T_SI2E , 0.0, 1.0); 
GDDCLI_Data(SPT_ADCSnapShot, ADC_SNAPSH0T_SIZE, adc_snap> ; 

GDDCLI_SendPacket (SPT_ADC_OOR, ( FLOAT3 2 ) ADC_GetOutOf RangeCount ( ) ) ; 

time = get_snapshot (ir, rd, raw_snapshot_size) ; 

/* validate this snapshot of data */ 

for (i = 0; i < raw_snapshot_size; i++) { 

if (ir[i] > 2.0) assert (0) ; 

if (ir[i] < -2.0) assert(O); 

if (rd[i] > 2.0) assert (0) ; 
^ if (rd[i] <-2.0) assert(O); 

store_raw(alg [0] .vp, ir, rd, raw_snapshot_size, time) ; 

compute_log(ir, log_ir, raw_snapshot_size) ; 
compute_log(rd, log_rd, raw_snapshot_size) ; 

for (i = 0; i < num_algorithms ; i++) { 

/* setup pointers for current algorithm */ 

vp = alg [i] .vp; 
vr = alg[i] .vr; 
vhp = alg [i] .vhp; 
vhr = alg[i] .vhr; 

/* allocate variable sized results arrays */ 

alg [i] . power = get_mem(vp->search_size * (vp->bin_count + 1) ) ; 
alg [i] .f spec = get_mem (vsi_spectral_size (vhp) ) ; 
alg [i] .pleth = get_mem (pleth_size (vp) ) ; 

power = alg[i] .power; 



fspec = alg [i] .fspec ; 
pleth = alg [i] .pleth; 

/* time stamp data */ 

vr->time = (FLOAT32) time / vhp->raw_sampling_rate; 

iawa_sat ( log_ir , log_rd, vp, vr) ; 

vsinc(ir, log_ir, log_rd, vp, vr, power); 
sat_f ilter (vp, vr) ; 

vsi_hr( log_ir, log_rd, vp, vr, vhp, vhr, fspec); */ 

vsi_hr (ir, log_ir, log_rd, vp, vr, vhp, vhr, fspec); 

hr_f ilter (vhp, vhr, vp, vr) ; 

vsinc_pl (ir, log_ir, log_rd, vr->art_omega_pr, vp, pleth); */ 
vsinc_j3l (ir, log_ir, log_rd, 
vr - >art_omega_j?r , 
vhp->hr_lambda, /* HR's lambda for plethy. */ 

vp, pleth) ; 

vr->Spa02_dual = 0.0; 

VSIDUAL_sat (vr->art_omega_pr, 
vp->dual_lambda, 
vp, 

ir, rd, log_ir, log_rd, 
&vr->Spa02_dual) ; 

GDDCLI_SendPacket (SPT_time_tick, (FLOAT32) time) ; 
GDDCLI_SendPacket (SPT_sat_f inal, vr->Spa02_f inal) ; 
GDDCLI_SendPacket(SPT_hr_final, vhr- >hr_f inal) ; 

} /* for */ 

store_j>arameters (num_algorithms) ; 
store_results (alg, num_algorithms ) ; 

/* send data */ 

f ree_mem ( ) ; 

mexit ( "single_step" ) ; 

static VOID 
test_pattern (void) { 

/* send a test pattern to verify graphing stuff */ 

} 

static VOID 
set_alg ( 

ALGORITHM_SPEC *alg, 

INT32 " num_alg) { 

SAT_PARAMS *vp; 
HR_PARAMS *vhp; 
CALIB_PARAMS * cp ; 

INT32 i; 

BUG1 (alg) ; BUG1 (num_alg) ,- 

for (i = 0; i < num_alg; i++) { 



vp = alg [i] . vp; 

vp->cov_limit 

vp- >num_points 

vp->anc_type 

vp->lambda 

vp->min_error 

vp- >anc_settle_time 

vp- >num_cells 

vp- >peak_limitl 

vp- >peak_epsilon 

vp- >:b P_ coe f 

vp->bp_size 

vp- >search_f irst 

vp- >search_step 

vp - > sear ch_s i z e 

vp- >bin_count 

vp->DVHS_init 

vp->DVHSO_init 

vp- >DVAST_timer 

vp - >DVAST_de 1 1 a 

vp - >DWST_t imer 

vp->DWST_delta 

vp- >max_sat_delta 

vp->sat_eqn 

vhp = alg[i] .vhp; 

vhp- >resolution 

vhp- >raw_sampling_rate 

vhp->min_hr 

vhp->max_hr 

vhp->nst 

vhp - >max_hr_de 1 ta 
vhp->hr enable 



= param_val [i] 
= (INT32) param_ 
= {INT32 ) paratn 
= param_val [i] 
= param_val [i] 
= (INT32)param_ 
= (INT32)param 
= param_val [i] 
= param_val [i] 
= post_coef; 
= post_size; 
= param_val [i] 
= param_val [i] 
= (INT32)param 
= (INT32)param 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= param_val [i] 
= (INT32)param 



[P_cov_limit] ; 
_val [i] [P_num_jpoints] ; 
_val [i] [P_anc_type] ; 
[P_lambda] ; 
[P_min_error] ; 
val [i] [P_anc_settle_tirae] ; 
_val [i] [P_num_cells] ; 
[P_peak_limitl] ; 
[P_peak_epsilon] ; 



[P_search_f irst] ; 
[P_search_step] ; 
_val [i] [P_search_size] ; 
_val [i] [P_bin_count] ; 
[P_dvhs_init] ; 
[P_dvhsO__init] ; 
[P_DVAST_timer] ; 
[P_DVAST_delta] ; 
[P_DWST_timer] ; 
[P_DWST_delta] ; 
[P_max_sat_delta] ,- 
i_val [i] [P_sat_eqn] ; 



1.0; 

(INT32)param_val [i] [P_sampling_rate] ; 
25.0; /* bpm */ 

225.0; /* bpm */ 
0.04; /* thresshold */ 

parara_val [i] [P_max_hr_delta] ,- 
(INT32)paratn_val [i] [P_hr_enable] ; 



vhp->hr lambda 




param_val [i] [P_hr_lambda] ; 


vp- 


>dual_lambda 




param_val [i] [P_dual_lambda] ; 


vp- 


>correlation 




param_val [i] [P_correlation] ; 


cp 


= alg'[i] .cp; 








cp- 


>rd_min_mA 




20.0; /* mA 


*/ 


cp- 


>rd_max_mA 




50.0; /* mA 


*/ 


cp- 


>ir_min_mA 




5.0; /* mA 


*/ 


cp- 


>ir_max_mA 




50.0; /* mA 


*/ 


cp- 


>target 




param_val [i] [P_cp_target] ; 


cp- 


>error_band 




param val [i] [P cp error band] ; 


cp- 


>ir_drive_l 




5.0; /* mA 


*/ 


cp- 


>rd_drive_l 




20.0; /* mA 


*/ 


cp- 


>ir_drive_2 




10.0; /* mA 


*/ 


cp- 


>rd_drive_2 




25.0; /* mA 


*/ 


cp- 


>settle_time 




(INT32)param val [i] [P 


_cp_settle_ 


cp- 


>ir_gain 




param_val [i] [P_cp_ir_gain] ; 


cp- 


>rd_gain 




param_val [i] [P_cp_rd_gain] ; 


cp- 






0.1; 





#define DATA_STEP 57 /* amount to scroll raw data */ 

VOID 

user_interf ace ( 

INT32 num_algorithms) { 

ALGOR I THM_S PEC *alg; 



INT32 



alg_num; 



/* current algorithm 



■I 





INT32 


param_num ; 


/* current parameter number 


*/ i 


615 


BOOL 


quit ■ 


/* quit program 






BOOL 


automation; 


/* quit automated session? 


'l 




INT32 


lsl_cmd; 


/* Low Speed Link command 


*/ ; 




FLOAT32 


*ir; 


/* for data preview 


*/ 


620 


FLOAT3 2 












*review size 
p — 


650 • 






INT 3 2 


ss size- 


/* sn h 

aps o size 




625 


INT32 










INT32 


j ; 








BOOL 


cal_ok; 








INT32 


time ; 








FLOAT32 


f ■ 






630 












BUG1 (nun 


_algorithms) ; 








assert (r 


um_algorithms <= 


MAX_ALGORITHMS) ; 





mark_mem ( ) ,- 

alg = get_mem(num_algorithms * sizeof (ALGORITHM_SPEC) ) ,- 

/* now allocate structures for each algorithm */ 

for (i = 0; i < num_algorithms ; i++) { 

alg[i].vp = get_mem (sizeof (SATJPARAMS) ) ; 

alg [i] . vr = get_mem (sizeof ( SATJRESULTS ) ) ; 

alg[i].vhp = get_mem (sizeof (HR_PARAMS) ) ,- 

alg[i].vhr = get_mem (sizeof (HR_RESULTS) ) ; 

alg[i].cp = get_mem ( sizeof (CALIB_PARAMS) ) ; 
^ alg [i] . cr = get_mem (sizeof ( CALIB_RESULTS ) ) ; 

/* init param structures */ 

for (i = 0; i < num_algorithms; i++) { 
for (j = 0; j < NUM_PARAMS; j++) { 
^ param_val [i] [j] = params [ j ] . def_val ; 

param_num = 0; 
alg_num = 0; 
quit = FALSE; 

ir = get_mem(preview_size) ; 
rd = get_mem(preview_size) ; 

new_file (0) ; 

set_alg(alg, num_algorithms) ; 
new_jpatient (alg, num_algorithms ) ; 

store ^parameters (num_algorithms) ; /* display initial parameters */ 
GDDCLI_FlushBuf f er ( ) ; 

CALIB_Init (alg [0] .cp, alg[0] .cr) ; 
do { 

GDDCLI_SendPacket (SPT_CalibrateMsg, 1.0); 
GDDCLI_FlushBuf f er ( ) ; 

cal_ok = CALIB_TwoPointCal (alg[0] .cp, alg[0].cr); 
GDDCLI_Data (SPT_CalibResults , CALIB_RESULTS_SIZE, 

(FLOAT32 *) (alg[0] .cr) ) 

} while (cal_ok == FALSE) ; 
GDDCLI_SendPacket (SPT_CalibrateMsg, 0.0); 



GDDCLI_FlushBuf f er ( ) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time); 
685 GDDCLI_FlushBuffer () ; 

presentation_f lag = FALSE; 

while (quit == FALSE) { 

690 

GDDCLI_HighLight (SPT_SystemParam, param_num, alg_num, 

param_val [alg_num] [param_num] ) ; 

GDDCLI_FlushBuf f er ( ) ; 
lsl_cmd = lsl_read(); 
695 GDDCLI_Normal (SPT_SystemParam, param_num, alg_num, 

param_val [alg_num] [param_num] ) ; 

GDDCLI_FlushBuffer() ; 

switch (lsl_cmd) { 
700 case NLT: /* left arrow */ 

if (alg_num > 0) { 
alg_num -= 1; 

break; 

705 

case NRT: /* right arrow */ 

if (alg_num < num_algorithms-l) { 
alg_num += 1; 

} 

710 break; 

case NUP: /* up arrow */ 
if (param_num > 0) { 
param_num -= 1; 

715 } 

break; 

case NDN: /* down arrow */ 

if (param_num < NUM_PARAMS - 1 ) { 
720 param nura += 1; 

} 

break; 
case NINS : 

725 f = params [param_num] . inc_val * 100.0; /* HGK */ 

goto kplus; 

case NPLUS: 

f = params [param_num] . inc_val ; /* HGK */ 

73 0 #if (PLATFORM == 1) 

if (shift_key() ) f*= 10.0; 

#endif 
kplus : 

735 param_val [alg_num] [param_num] += f; 

if (param_val [alg_num] [param_num] > params [param_num] .max_val) 
param_val [alg_num] [parara_num] = params tparam_num] . max_val ; 
break; 

740 case NDEL : 

f = params [param_num] . inc_val * 100.0; /* HGK */ 

goto kminus; 

case NMINUS: 

745 f = params [param_num] . inc_val ; /* HGK */ 

#if (PLATFORM == 1) 

if (shift_key() ) f*= 10.0; 

#endif 



param_val [alg_num] [parara_num] -= f; 

if (param_val [alg_num] [param_num] < params [param_num] .min_val) 
param_val [alg_num] [param_num] = params [param_num] .min_val; 
break; 



case NRTUP: /* page up */ 

move_snapshot ( -DATA_STEP) ; 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg[0] .vp,ir, rd, preview_size, time); 

GDDCLI_FlushBuf f er ( ) ; 

break; 

case NRTDN: /* page down */ 

move_snapshot ( DATA_STEP) ,- 

time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp,ir, rd, preview_slze, time) ; 

GDDCLI_FlushBuf f er ( ) ; 

break; 

case SCR: /* continuous operation */ 
set_alg(alg, num_algorithms) ,- 

ss_size = compute_snapshot_size (alg, num_algorithms) ; 
'* calibrate2 (system_params . cal_target) ; */ 
reset_file_num() ; 
new_file(0) ; 

automation = TRUE; /* "continue automated session' */ 

do 

{ 

time=0; 

new_patient (alg, num_algorithms) ; 

do 

{ 

single_step(alg, num_algorithms, ss_size) ; 
move_snapshot (USR1_AUT0_STEP) ; 

/* mating get_snapshot ( ) is in single_step . . */ 
if (lsl_status() ) 

automation = FALSE; /* discontinue session */ 

while (automation && ! end_of _data ( ) ) ; 
if (automation) 

GDDCLI NewFileO ; 

} 

while (automation && ! end_of _f iles ( ) ) ,- 
break; 

case SCS: /* single step */ 

set_alg(alg, num_algorithms) ; 

ss_size = compute_snapshot_size (alg, num_algorithms) ; 
single_step (alg, num_algorithms, ss_size) ; 
break ; 

case SCC: /* "C" calibrate command */ 



case SCT: /* "T" test pattern command */ 
test_pattern() ; 



case NESC: /* escape */ 
quit = TRUE; 
break; 

case SCF: /* "F" new file command */ 
new_f ile (1) ; 



time = get_snapshot (ir, rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time); 
GDDCLI_FlushBuf f er ( ) ; ; 
break; 

case SCB: /* "B" back file */ 
new_f ile (-1) ,- 

time = get_snapshot ( ir , rd, preview_size) ; 

store_raw(alg [0] .vp, ir, rd, preview_size, time) ; 
GDDCLI_FlushBuffer() ,- 
break ; 

case SCH : 

param_val [alg_num] [P_num_cells ] =2 ; 
param_val [alg_num] [P_search_size] =25; 



#if (PLATFORM == 1) 
case SCW 



default: /* ignore unknown commands */ 
break; 



} /* switch */ 
} /* while */ 



f ree_mem ( ) ; 
} /* user_interface */ 



